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Abstract. Basing on the phenomenological model of the microwave response of thin 
superconducting films with trapped magnetic flux elaborated in Luzhbin et al 2005 Supercond. 
Sci. Technol. 18 1112-7, detailed analysis of the ac linear response of a superconducting 
strip-line resonator cooled down in FC and ZFC regimes is presented. Special emphasis on 
a nonuniform distribution of the trapped vortices and on the magnetic field and temperature 
dependence of the phenomenological model parameters such as critical current density, vortex 
pinning and viscous motion coefficients is made. For a ZFC case realistic trapped magnetic 
field distribution based on the critical-state model for a vortex penetration is used. Within the 
framework of the used model the unusual multi-peak dependence of the surface resistance of 
HTS films observed experimentally is explained as well as clarified conditions for its appearance. 



1. Introduction 

The phenomenological theory of the linear response of type-II superconductors (SC) in the mixed 
state developed in the last decades El El El El describes the experimental data sufficiently well. 
Being based mainly on the idea of the induced current and its effect on the quasi-particles due 
to vortex motion, it can be formulated in terms of a change of the local effective resistivity of 
the superconducting material. Introducing into the model description various phenomenological 
parameters (e.g. pinning constant a^, viscosity r\, creep rates r, tq [21 Ej, critical current 
density J c , elastic moduli of the flux line lattice Cu, i = 1,4, 6, etc) allows us to simultaneously 
incorporate into the theory various effects of the external magnetic field and the temperature 
starting from the ordinary flux creep to the phase transitions in the vortex ensemble. 

However, a number of experiments shows unusual quasi-periodic dependence of the surface 
impedance of HTS films on the temperature or the external dc magnetic field (El 13 El) which 
cannot be unambiguously explained within the frameworks of the existing theoretical models. 

In the current report the theoretical model of the linear response of HTS thin films with 
trapped vortices developed in [2] is used to model the dc magnetic field and the temperature 
dependence of energy losses Pdi s (B,T) connected with the oscillation of Abrikosov vortices in 
thin SC films cooled down in (FC) and without (ZFC) external dc magnetic field. It is shown 
that, under the general assumptions on the magnetic field and temperature behavior of J c , ax 
and r\, the quasi-periodic form of -Pdis(^) dependence neatly corresponds to the ones obtained in 
00 El- By varying the shape of J C (B,T), a^(i?,T) and rj(B,T) dependences, the conditions 
for the appearance of nonmonotonic parts on Pdi s (B,T) curve is also clarified. 



2. Theoretical consideration 

Let us consider a thin SC film (with thickness d < X and width 2w S> 2X 2 /d) with a 
trapped magnetic flux. It is supposed that all physical quantities (current, magnetic field, 
etc) are distributed along the y-axis only and averaged over the film thickness, i.e. J(y) = 



The trapped magnetic flux is produced by Abrikosov vortices distributed along the film width 
with some magnetic field profile B(y), which is proportional to the local density of vortices. The 
profile B(y) in turn depends on the cooling prehistory of the film and can in principle be found 
by means of magneto-optical, scanning Hall probe microscopy, or Lorentz microscopy methods. 

Upon the application of an external ac current J r f (y) ex.p(—iu)t), vortices are forced to oscillate 
around the equilibrium positions with the same frequency uj. The vortex displacement field u(y) 
can be found from the force balance condition [21 El, which in our case has the form 



where a = — ilo/loq), ujq = oil/t] is the dep inning frequency that lies in the range 

10 10 — 10 12 s _1 JS], and $o is the A ux quantum. Equation (J^) models the well known ac 
current concentration in superconducting strip-line resonators at the film edges by a function 
with the square root singularities at the edges of the film, which is a rather good approximation 
jlUllllj . Moreover, since equation (^Q) has the finite solution over the whole interval —w < y < w, 
it proposes a very realistic model of the linear response of thin SC films with trapped vortices. 

Having found the vortex displacement field u(y) as the solution to one can calculate the 
response of the film to the ac current; in particular, the ohmic losses Pdis per unit length of the 
film during the period of the ac field oscillations can be found according to the definition ^2] 



where all time dependent variables are supposed to have the form exp(— iut), g is defined by the 
geometry and operation regime of the resonator, y e dge = w — 2A 2 /d, and P e dgc corresponds to 
the ohmic losses at the edges of the film, i.e. for y e dge < \y\ < w i J{u) 1S the total current in the 
resonator, and E{y) is the ac electrical field generated by the current J(y). In (j2j) the fact that 
the real current distribution in the resonator is obtained by the current cut off at a distance of 
order 2X 2 /d from the edges is taken into account. 

The general idea of the numerical solution of equation Q that gives the approximate solution 
with an arbitrary specified precision can be summarized as follows |13j . 

(i) For any smooth function B(y) that has no zeros in the interval [—w, w], the function f{y) = 
B{y)u{y) can be well approximated by a polynomial of degree N, given by f(y) ~ J2jLo a jV 3 ' ■> 
where a,j are the unknown complex coefficients, y = y/w, and f'(y) — J2jLoj a jy^~ 1 - 

(ii) Let us write B{y) = B dc f B (y) and J r f(y) = J r f,o/j(y), where B dc and J rf)0 are some 
constants that have dimensions of magnetic field and current density, respectively, and the 
functions fs and fj define the dc magnetic field profile and the ac current distribution over 
the film width, respectively. In general, /s(y) also depends on the external dc field |14j . 
but in the simplest case of the uniform vortex distribution fg is a function of y only, and 
B dc h as the meaning of the external dc field (e.g. for FC case). For the linear response it is 
always possible to write a,j = dw^oJ^^Cj, where Cj are dimensionless unknown coefficients. 
Multiplying both sides of equation Q by B(y) yields the equivalent equation 



L d/2 J{y,z)dz/d. 
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where S(u>) = dw/ \q(uj), Xq(uj) = B& c <&q/ hqo. is the Campbell penetration depth [2j EJ, 
and Vj (y) = f\ (^VR 1 ^ /(C - V). 
(iii) Taking the values of all functions in © at the basic points yo, y\, ... y^ which are the 
zeros of the Chebyshev polynomial Tv + i(y) = cos((iV-|- 1) arccos(y)), one obtains a system 
of linear equations for determining the unknown coefficients Cj . 

Since the parameters defining the response of resonators in this approximation (i.e. fj(y) 
and are supposed to be known from elsewhere, this model is directly applicable to any 

type of thin film resonator. The quasi-periodic form of Pdi s (S(u>)) dependence that follows from 
equations (|T|). © (see figure^ A)) is a manifestation of the dimensional effect associated with 
quantization of the displacement field u(y) on the width of the film 

To model the temperature or the dc magnetic field dependence of the microwave response 
in the framework of this model, one should take into account the following factors: the 
temperature/magnetic field dependence of the Jb profile (e.g. due to flux creep), and the 
dc magnetic field and temperature dependence of the phenomenological parameters ocl, f] and 
J c , which is supposed to be known from elsewhere, see, e.g. [TBI HEIIT7] . 

For the sake of an example let us consider the temperature dependence of the energy losses 
for a film cooled down in FC and ZFC regimes in an uniform external magnetic field B^ c . For 
the FC case the uniform distribution of the trapped vortices is a rather good fit, i.e. B(y) = B^ c , 
whereas for the ZFC case the critical state model ^1] is the best approximation. The original 
critical state model is modified in the following way: B(y) = B rcmn for |y| < b, and \ J{y)\ = J c 
for b < \y\ < w with b = wj cosh(7ri? ( j c //io J c ) and B rcmn "C B& c being some remnant field. This 
modification allows us to use the foregoing approach to model the ac response of the film in 
ZFC state and, as long as B remn <C -Bdo it does not affect the results; on the other hand, such 
modification is justified by the fact that in a lot of experiments some unscreened magnetic field 
of order of several oersted is always present jH]. In both cases, the following temperature ^3] 
and the dc magnetic field JS] dependence of oti and rj is supposed: 

a L (T) = «l(0)(1 - tfl\l + tf exp(-T/T ) , a L (B) ex B% , 
rj(T) = 77(0X1 - i 2 )/(l + t 2 ) , r,(B) oc 2?£ , 

where t = T/T c , T c = 88K, typical values of a L (0) ~ 3 x 10 5 N m~ 2 and r?(0) ~ (0.2 - 1.2) x 
1CP 6 N s m -2 , To is a characteristic fluctuation temperature [^I15j (with To ~ 28K for YBCO), 
n a ~ 0.55, n r] ~ 1.85 16 j. Note that the given values of O£,(0), r/(0) as well as their temperature 
and field dependences (@J are the best fit for the numerous experimental data obtained in a 
number of experiments on YBCO |15[ I16j. whereas theoretical predictions vary from the field- 
independent rj for the isotropic s-wave pairing to n v = 1/2 for the d-wave pairing, and from the 
field-independent oil for the individual pinning of vortices to n a = 1 for the collective pinning. 

The temperature and the magnetic field dependence of J c is well described by the following 
empirical formula [H|: J C (B,T) = J c0 (B)(l-t) £ , where J c0 (B) = J c0 (0)/3 ln(B*/B) in the range 
of interest, B* ~ -B c ff(l — t). J c o(0), e, (3 and B e g are some empirical parameters, which can be 
independently obtained from transport, ac susceptibility and magnetization measurements. For 
YBCO films, the following values are typical: (5, e are of order of unity, B e g of order of several 
tesla, and J c o(0) of order of several MA cm -2 [T7j . 

3. Discussion 

The energy loss dependence on the temperature Tdis(T) for a film in FC and ZFC states is shown 
in figure^B). All the parameters ol, r], and J c are supposed to be temperature and magnetic 
field dependent according to ((H) with ol(0) = 10 5 N m~ 2 and 77(0) = 0.2 x 10~ 6 N s m~ 2 . 

It is found that the exact form of Pdi s (B, T) curve and the peak positions are mainly defined 
by «l(0) and J c o(0) as well as the exact form of J C (B,T) dependence, while the quantities rj 




Figure 1. Vortex contribution to the ohmic losses in the film (A) as a function of the parameter 

So = dw/\^j{u> = 0), Q = lo/uq, and (B) as a function of the temperature for FC ( ) and 

ZFC ( ) regimes for the frequency lo/2it = 1 GHz. To avoid crowding, ZFC curves are 

multiplied by factor 4. 



and too define the peak amplitudes. Decrease of ol(0) and J c o(0) shifts the position of the peaks 
to the lower temperatures and decreases their amplitude, so that the peaks cannot be detected. 

Other crucial factors are the geometrical parameters of a resonator, i.e. its thickness d and 
width 2w. As it is seen from figure ^A), the most pronounced peaks in the curve -Pdis(£b) are 
located in the region of tens of So; thus, only if the dimensions of a resonator lie in the range 
corresponding to 5 < So < (20 — 30) for the rather low frequency, the multi-peak picture similar 
to the ones obtained in jSJ[71|S] can be experimentally observed. 

In that way, the possibility to experimentally observe the quasi-periodic (multi-peak) 
dependence Pdi s (T) crucially depends on the transport properties of HTS films, which can 
be characterized by J c and ol and in turn are determined by the conditions and method of 
preparation of the films, as well as on the geometrical dimensions of SC resonators. 
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